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The cusp effect mK ^3>n and data on Kg/[ decays allow one to extract experimental information 
on the elastic nn scattering amplitude near threshold, and to confront the outcome of the analyses 
with predictions made in the framework of QCD. In my talk I pointed out that these predictions 
concern an isospin symmetric world, while experiments are carried out in the real world, where 
isospin breaking effects - generated by electromagnetic interactions and by the mass difference 
of the up and down quarks - are always present. I discussed the corrections required to account 
for these, so that a meaningful comparison with the predictions becomes possible. In particular, 
I pointed out that there is a spectacular isospin breaking effect in KeA decays, overlooked so far. 
Once it is taken into account, the previous discrepancy between NA48/2 data on KeA decays and 
the prediction of nn scattering lengths disappears. 
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1. Introduction 

Chiral perturbation theory (ChPT) [|l|], combined with Roy equations, allows one to make 
very precise predictions for the values of the threshold parameters in elastic nn scattering - see 
Colangelo's contribution to these proceedings for a status report [||]. Several experiments allow one 
to confront these predictions with the data: i) n^n^e^Ve decays [||, ^, ||], ii) the pionium 

lifetime, measured by the DIRAC collaboration [^, and iii) the cusp effect in A' ^ 371 decays, 
investigated by the NA48/2 collaboration 

The experiments performed by the NA48/2 collaboration have generated an impressive data 
basis, as a result of which the matrix elements of Ke4 and K ^3n decays can be determined with 
an unprecedented accuracy. The interpretation of these measurements was the topic of my talk. 
In particular, I pointed out that the theoretical predictions and the measurements are performed in 
two different settings: the predictions concern pure QCD, in the isospin symmetry limit m„ = m^, 
with photons absent. To be more precise, the convention is to choose the quark masses and the 
renormalization group invariant scale of QCD such that the pion and the kaon masses coincide 
with the values of the charged ones, and the pion decay constant is Fji = 92.4 MeV. [I do not 
specify the masses of the heavy quarks, because in the present context, their precise values do not 
matter.] I refer to this framework as a paradise world. 

On the other hand, experiments are all carried out in the presence of isospin breaking effects, 
generated by real and virtual photons, and by the mass difference of the up and down quarks: this 
is the real world, described by the Standard Model. We are thus faced with the problem to find the 
relation between quantities measured in the real world, where isospin breaking effects are always 
present, and the predictions made in the paradise world. 

Remarks: The situation is quite involved: 1. The above mentioned three processes do not take 
place in the paradise world: the kaon is stable, whereas pionium cannot form, because electromag- 
netic interactions are absent. 2. In the real world, the processes i) and iii) do not occur either, 
because decays that involve external charged particles have zero probability to happen without the 
emission of photons. 3. Lattice calculations are not confronted with this problem: the evaluation 
of e.g. Tin or nK scattering lengths ^ can be performed directly in the framework of pure QCD, 
at mu = md: isospin breaking effects can be excluded by fiat. I refer the interested reader to Colan- 
gelo's contribution [^] , and to Ref. [Q] for a short review of lattice results concerning nn scattering 
lengths. 



2. The cusp inK -^3n 



Cabibbo and Isidori pointed out | [To| ] that cusp effects in the decay n^n^n^ may allow 

one to measure 5-wave tttt scattering lengths, and thus to check the precise predictions made in 



the framework of chiral perturbation theory, see also [11]. A first data analysis based on this 



observation has appeared [12] already some time ago, and the latest developments were discussed 
at this conference [^. 

We have started to investigate this process in the framework of non-relativistic quantum field 



theory in Ref. [13]. The status of our calculation is as follows. 
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Our amplitude agrees with the one of [10, 11] near the cusp, whereas it differs away from 



it. The NA48/2 collaboration at CERN has analyzed data with our amplitude as well. The 



outcome for the scattering length combination oq — a2 is in agreement [ ]14[ ] with the use of 



the amplitude from ||10|]. 



We have evaluated the decay amplitude for the neutral kaon, n^n^n^,n^n'^n^, have 

written a pertinent FORTRAN code for the decay probability and made it available to the 
NA48/2 collaboration |[T^]. In the channel with 3 neutral pions, a cusp occurs as well. In 
this case, its strength is suppressed [Q]. Nevertheless, the huge amount of data available at 
NA48/2 should make it possible to extract the same combination of scattering lengths as 
well, although with less precision [^. 

It remains to calculate electromagnetic corrections to these decays. Because we make use 
of a quantum field theory framework, our method allows us to perform these calculations 
in a straightforward, yet tedious manner. We have investigated the role of the formation 
of pionium near the cusp and have identified an interval around the cusp, where radiative 



corrections can be calculated reliably, in a systematic manner [15]. [In Ref . [16], the structure 
of the cusp in the presence of re-summed Coulomb ladders is investigated as well. In case 
that the same approximation is made, the so induced modifications of the decay distribution 
agree in the two frameworks. Furthermore, pionium production was calculated in [|l^], and 
compared with data by the NA48/2 collaboration 



Presently, we investigate the emission of real photons []15|]. We evaluate the infrared diver- 
gent contributions in ^ —> 371(7) in dimensional regularization, and organize the remaining 
integrations over phase space in such a manner that the calculations can be performed numer- 
ically reasonably fast, so that a corresponding FORTRAN code can be used by the NA48/2 
collaboration for performing fits. 



3. decays 
3.1 General 

In the NA48/2 experiment, the general purpose Monte Carlo program PHOTOS [jl^] is used 
to calculate electromagnetic corrections. In addition, the Sommerfeld factor is applied, to account 
for the Coulomb interaction between charged particles [|l9|]. 

In my talk, I pointed out that in these prescriptions to perform radiative corrections, one spe- 
cific mechanism is not included. Namely, the kaon may decay first into a neutral pion pair, that then 
annihilates into two charged pions, or first decay into a charged pion pair, that then re-scatters. In 
the real world, the neutral pion mass is smaller than the charged one by about 4.6 MeV^ As a result 
of this, the two contributions to the decay matrix element have a different holomorphic structure: 
the neutral (charged) pion loop generates a branch point at Sji = 4M^o (at s = 4M^), and the phase 
of the relevant form factor is affected with a cusp, and does not vanish at the threshold s = 4M^. 
In my talk, I discussed the case of the scalar form factor {n^ \uu + dd\0) , where the same effect 

' Throughout the text, I use the symbols Mjc (M^jo) for the charged (neutral) pion mass. 
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occurs. Assuming universality, I applied the result to Ke4 decays [20|. Meanwhile, we have worked 
out [ 21 ] the effect for Ke4 decays as well. I describe the outcome of the calculation shortly here, 
and refer the interested reader to our forthcoming publication for details [^. 

3.2 Partial wave expansion: isospin symmetry limit 

The matrix element for (/?) {pi)n^ {p2)e^ {pe)Ve{pv) is 



V2 



v:Mpv)y''{i-Y5)v{pe){v^-A^), 



where the last factor denotes hadronic matrix elements of the strangeness changing (vector and 
axial vector) currents, 

y^i-A^ = (7r+(pi)7r"(7?2)out I {sy^u - sy^jsu) \ K+{p)). 

In the following, I concentrate on the matrix element of the axial vector current, because it carries 
information on the nn final state interactions and, in particular, on the Tin phases. One decomposes 
into Lorentz scalars, 

Ali = -^Tr [{Pi+ P2)tiF + {P\- P2)tiG + {Pe + Pv)^K] . 

Mk 

The form factors F,G,K are, holomorphic functions of the three variables 

% = {Pl+P2f,t = {pi -pf,U= {p2-pf ■ 

Sometimes, it is useful to use instead 

Sn = {P\ + Plf , Si = {pe + Pvf , COS , 

where dj^ is the angle of the 71+ in the CM system of the two charged pions, with respect to 
the dipion line of flight in the rest system of the kaon [ p2[ ||]. In the isospin symmetry limit, 
one identifies the nn phases in the matrix element in a standard manner, by performing a partial 
wave expansion, and using unitarity and analyticity, although, in the present case, this is a slightly 
intricate endeavor [^]. It is useful to introduce a particular combination of form factors, 

{Ml-s^-se)o 
^1 = ^ + T7772 ^^^^^^G. 



Here, a = a/1 — AM'^/sj^, and X{x,y,z) is the triangle function. The form factor F\ has a simple 
partial wave expansion, 

k>l 

In the low-energy region Sj^ < 16M^, fk carry the nn phases |Q in the pertinent isospin channel. 
In the following, I consider the lowest partial wave f{sji,S(). In the interval 4M^ <Sji< 16M^, its 
phase coincides with the isospin zero S'-wave phase 5q in elastic nn scattering, 

f+ = e''^'f^J±=f{s^±ie,s,). (3.1) 
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Figure 1: Some of the graphs that contribute to the matrix element of the axial current at tree and one- 
loop order The filled vertex indicates that the axial current also couples to a single kaon line. That graph 
contributes to the form factor K. There are many additional graphs at one-loop order, not displayed in the 
figure. 

It is instructive to calculate the form factors in chiral perturbation theory and to verify that F\ indeed 
has the behaviour just discussed. For this, it is sufficient to consider the effective Lagrangian 

where the covariant derivative D^U contains the external vector and axial vector currents, and 
^ = dmg{m,m,ms). Some of the graphs that contribute at tree level and at one loop are displayed 
in Figure I. The result is [||] 



with 



A(5;,) = -}-^{2s^-Ml)J{s„) 



l6n^J{s^) = o[\a ^-p^ + inj +2, s„> AmI . 

Here, Mji (Fq) denotes the pion mass (pion decay constant), at leading order in the chiral expansion. 
The quantity H{sjc,se) is real in the interval of elastic nn scattering. It is now seen that / indeed 
has the property ( [3.1[ ) at this order in the low-energy expansion, with 

This is the phase of the isospin zero 5-wave, in tree approximation. 

3.3 Partial wave expansion: the real world 

In reality, experiments are not carried out in the paradise world just discussed: we have not 
included so far photons, nor did we consider isospin breaking effects generated by different up and 
down quark masses. Here, I investigate these effects in several steps: 
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Figure 2: The contribution from 7Z ~rj mixing, at leading order in — mi,. The filled square denotes the 
vertex from tt" — T] mixing. Otherwise, the notation is the same as in Figure 1 . 



i) I assume that the manner in which real and virtual photons are treated in the analysis of the 
NA48/2 experiment (PHOTOS + Sommerfeld factor) is a decent approximation to the effects 
generated by soft photons. 

ii) This procedure misses the effects generated by the pion and kaon mass differences, and by 
the quark mass difference — m„. These must therefore be taken into account separately. 

iii) ChPT is the appropriate tool to evaluate these contributions. 

iv) I assume that PHOTOS+Sommerfeld factor + mass effects provide a good approximation to 
the full isospin breaking contributions. 



Remark: One may envisage a more ambitious procedure [25], by working out the relevant matrix 



elements in the framework of ChPT including photons and leptons [26], and then constructing a 



new event generator, to be used in the analysis of decays. [A one-loop calculation was already 



performed in [27]. It needs to be checked, and brought into a form which is suitable for the present 
purpose.] Eventually, such an analysis might lead to an improved algorithm, but I consider it a long 
term project. End of remark. 

According to iii), we simply need to perform a ChPT calculation of the effects generated by 
the mass differences. This is rather easy to achieve at one-loop order: one adapts the quark mass 



matrix, diag (m„,md,m.v), and enlarges the Lagrangian ^2 [28], 



if2 ^ ^2+C{QUQW) , Q = ^diag(2, -1, -1) , (3.3) 
where C is a new low-energy constant, that breaks the isospin symmetry of the meson masses: 



The effect of the replacement ( |3.3| ) is twofold: first, as just mentioned, the meson masses split. 
As a result of this, the loop contributions in Fig. [l]b),c) have a different threshold, and the phase of 
the form factor / generates a cusp. Second, in addition to the graphs displayed in Figure |l], there 
is a new contribution shown in Figure ^: the kaon interacts with the axial current to generate a 
71° T] intermediate state. Because m„ 7^ m^, the T7 can transform back into a neutral pion, that then 
re-scatters with the second neutral pion into a charged pion pair. 



Working out the relevant diagrams, one finds that the phase ( |3.2| ) becomes in the elastic region 



5,? - 5 = I (4A, + s,)a + - m2„) (1 + aoj , (3.4) 
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with 



m 



md - niu 



The one-loop expressions for the form factors F, G given in Refs. [ ]27| ] contain the effects consid- 
ered here, up to terms of order (XQEoi^d — niu)- 
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Figure 3: The isospin breaking correction that must be subtracted from the phase 5 measured in K^n decays. 
The width of the band reflects the uncertainty in the ratio R. 



I consider the result ( |3.4| ) to be very interesting, for the following reasons. First, due to the 
presence of the phase space factor Oq, the phase 8 does not vanish at the threshold Sji = 4M^. This 
unexpected behaviour of the phase is the cusp effect already experienced in ^ ^ 371 decays, with 
the role of neutral and charged pions interchanged. Second, the difference 5 — 5q is positive for 
above the threshold, and even increases at large Sj^, 



3% 1 



6471^(2 R 



+ 1 



> 1 . 



We now come to the main point. According to point iv) above, it is the phase 8 that is measured 
in KeA decays (up to contributions from higher orders in the chiral expansion). Therefore, before 
comparing the phase so determined with ChPT predictions, one has to subtract from the measured 
phase the (positive) difference 8 — 8^, because 8^ = 8 — {8 — 5q ). In Figure ^ we display this dif- 
ference in the relevant decay region, for /? = 37 ±4^. The width of the band reflects the uncertainty 

^This value for R should be considered as preliminary ~ it was used in my talk for illustrative purposes. A more 
refined estimate will be provided in [^]. Of course, the conclusions to be drawn from the isospin breaking effects 
considered here will not change. 
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in R. [As I pointed out in my talk, two-loop contributions are modest in the analogous case of the 



scalar form factor of the pion []21[].] It is seen that the isospin correction is quite substantial - well 



above the uncertainties quoted for the measured phase [g]. [In Ref. []29|], the cusp in Ke4 decays 
was investigated as well. The expressions presented there do not agree with (3^), because these 
authors do not take into account derivative couplings of the nn amplitude, as is dictated by chiral 
symmetry.] 

Colangelo has performed fits to Ke4 data, with and without isospin breaking corrections ap- 



plied. It turns out that the former discrepancy [ |30| ] of the NA48/2 data with the prediction 



disappears, once isospin breaking effects are taken into account in the manner just described, see 
[^]. Colangelo also shows that the former agreement between the chiral prediction and the E865 
data [Q] becomes marginal. This is an issue that should be understood, because it is independent of 
the special effects considered here. On the other hand, since the NA48/2 data are so precise, they 
will dominate the E865 result, in any case. 

Note added: Isidori has meanwhile worked out soft photon corrections for — > tt+tt+tt" 



[31] 
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